Low-temperature quantum fluctuations in overdamped ratchets 
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At low temperatures and strong friction the time evolution of the density distribution in position 
follows a quantum Smoluchowski equation. Recently, also higher-order contributions of quantum 
fluctuations to drift and diffusion coefficients have been systematically derived. As a non-trivial 
situation to reveal the impact of subleading quantum corrections and to demonstrate convergence 
properties of the perturbation series, directed transport in ratchets is studied. It is shown that the 
perturbation series typically has a non-monotonous behavior. Depending on symmetry properties 
higher-order contributions may even compensate current reversals induced by leading quantum 
fluctuations. This analysis demonstrates how to consistently treat the dynamics of overdamped 
quantum systems at low temperatures also in numerical applications. 
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PACS numbers: 05.40.-a,05.60.Gg,05.70.Ln,73.23.-b 
I. INTRODUCTION 

Directed transport in ratchet-type of force fields has 
attracted substantial research for more than a decade 
now as one prototype of Brownian motor (see e.g. QUI])- 
These engines perform work by extracting energy from 
out-of-equilibrium fluctuations of the environment. Re- 
alizations can be found in a variety of contexts ranging 
from biological systems to designed mesoscopic devices 
such as electrical circuits containing Josephson junctions. 
The corresponding classical theory is well established 
and based on Langevin equations and, equivalently, on 
Fokker-Planck equations for the phase-space dynamics. 

Much less is known for the corresponding quantum me- 
chanical problem, mainly due to the fact that in this 
regime the interplay of quantum fluctuations in the sur- 
rounding heat bath and system dynamics generates long- 
range retardation effects in time so that a simple time 
evolution equation for the density distribution does not 
exist. Progress has been made particularly in the adia- 
batic limit for a dichotomous noisy drive and for tight- 
binding systems 0-11] • A considerable simplification oc- 
curs in the strong-friction range. Then, the time scale for 
relaxation in position by far exceeds all other time scales 
and a time-local evolution equation for the marginal po- 
sition distribution can be derived from the formally ex- 
act path-integral expression. This so-called quantum 
Smoluchowski equation (QSE) |9j applies also to the low- 
temperature regime, where the energy scale for friction 
fvy (friction constant 7) is much larger than the thermal 
energy scale kg,T. Conceptually, this domain is far from 
being classical which requires just the opposite condition 
jh <C &bT, even though quantum fluctuations can be 
treated perturbatively. Their leading contributions are 
of order 111(7) ll an d typically have strong impact in bar- 
rier escape problems. 



This is particularly true for transport in ratchets as 
shown in detail in a series of papers Net cur- 

rents are extremely sensitive to variations of the poten- 
tial topology so that weak quantum fluctuation may e.g. 
even generate current reversals. What has not been done 
though, is a systematic analysis of higher-order quantum 
corrections. This is the goal of the present work based 
on the general perturbative approach for drift and dif- 
fusion coefficients in the QSE recently developed in [13| . 
The naive expectation is that subleading quantum fluc- 
tuations lead to monotonous convergence properties. In 
contrast, our results indicate that typically this is not 
the case and even current reversals may be completely 
compensated by higher-order contributions in the per- 
turbation series. It turns out that a consistent treatment 
has to take into account orders in perturbation theory, 
which carry information about non-local characteristics 
of the ratchet potential. More generally, the results of 
this study show how subleading quantum fluctuations 
must be consistently incorporated in an extended QSE 
to describe the real-time dynamics of strongly condensed 
phase systems at very low temperatures. 



II. QUANTUM SMOLUCHOWSKI EQUATION 

In preceding publications 0, [l3j], we have considered 
quantum dissipative systems in the limit of strong fric- 
tion, in which a typical damping constant 7 by far ex- 
ceeds a characteristic frequency lu of the isolated sys- 
tem. In the corresponding classical domain, the so-called 
Smoluchowski regime is characterized by a separation of 
time scales, whereas the momentum decays on the time 
scale I/7, while its position relaxes on a much larger 
scale 7/w 2 [I3 - [l6j . The relevant dynamics of this latter 
process is described by the Smoluchowski equation (SE) 
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for the marginal position distribution P(q, t) of a 
particle with mass m in a potential V(q) surrounded by 
an environment at inverse temperature (3 = \/(k&T). 

Quantum mechanically the situation is much more in- 
volved. The intrinsic thermal time-scale for fluctuations 
of the bath Hf3 comes into play such that even an envi- 
ronment with white-noise characteristics in the classical 
domain exhibits colored-noise fluctuations at lower tem- 
peratures. Accordingly, simple dynamical equations for 
the reduced density matrix (results from the full density 
matrix of system+reservoir after tracing out environmen- 
tal degrees of freedom) do not exist. For strong friction, 
an effectively Markovian diffusion-equation for the diag- 
onal part of the density, the position distribution P(q, t), 
can thus only be derived if the time-scale separation is 
extended to include also the thermal time-scale. The cor- 
responding QSE has been obtained in leading order in 
the quantum fluctuations in [9| and since then has been 
applied to a variety of transport problems (see [l8| and 
references therein). Recently, subleading quantum cor- 
rections to this original QSE have been systematically 
derived within a type of semiclassical analysis in [131 ] to 
which we refer for further details. Here we only recall the 
main features. A general QSE can be cast in the form 
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with drift and diffusion coefficients D\ and D 2 , respec- 
tively. In principle, these coefficients must be derived 
from the full path-integral dynamics in real time. How- 
ever, in case of a system with time-independent poten- 
tial V(q), it turns out that they can be determined by 
thermal-equilibrium properties. The corresponding Eu- 
clidian path-integral approach shows that in the strong- 
friction regime the thermal distribution in position Pp (q) 
is obtained as 



P (q) = Z- 1 F(q)e 



(2) 



with a minimal Euclidian action ip(q), a fluctuation pre- 
factor F(q), and a normalization Z. Then, drift and 
diffusion coefficients in (JTJ can be read off from these 
quantities as 
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Note that this form of the coefficients guarantees that 
spurious equilibrium currents are avoided in each order 
of perturbation theory [l(], EH. The factor d(q) in ([1]) 
captures dynamical and additional thermal corrections 
and is obtained from the full time-evolution of the re- 
duced density matrix. It turns out that up to the order 
in perturbation theory, where analytical results are avail- 
able, dynamical corrections are completely classical and 



are at most of order uj 2 /j, while additional thermal fluc- 
tuations are fixed through the exactly solvable case of 
a harmonic system We note in passing that the 

current operator in ([T]) may also be cast in the standard 
Fokker-Planck form Di + d q D 2 with D 1 = dD 1 - (d q d)D 2 
and D 2 = dD 2 . 

As will be seen in detail below, quantum corrections 
in the coefficients of the QSE carry higher than first or- 
der derivatives of the potential V(q). This is attributed 
to the general tendency of quantum fluctuations for de- 
localization in position and in momentum in accordance 
with the uncertainty principle. For externally driven sys- 
tems with time-dependent potentials V(q, t) the semiclas- 
sical evaluation of the full path-integral dynamics reveals 
that the results of the static case are simply generalized 
by replacing V(q) — > V(q,t) if the external driving hap- 
pens to be sufficiently smooth and if its typical frequency 
O respects the time-scale separation I/7, 7i/3, ft/3 < l/fl 

mug. 

To analyze the impact of higher-order quantum cor- 
rections in more detail, we concentrate in the sequel 
on the low-temperature Smoluchowski regime given by 
1/7 <C ft/? -C A systematic perturbative treat- 

ment of the thermal distribution ^ is then based on a 
power series in the small parameter 
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where \& denotes the digamma function and Ce is the 
Euler-Mascheroni constant. In leading order (order A) 
the current operator is found to read 
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which gives the QSE already used previously, see e.g. 

m- 

In next order (order A 2 ), the path integral for the equi- 
librium distribution Pp(q) is solved by taking into ac- 
count local harmonic properties of the potential around 
the endpoint q of the closed minimal-action paths. 
Hence, we find 
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For even higher-order terms, local anharmonicities are 
considered as well. Analytical results are available up 
to contributions of order A 2 /L with an anharmonicity 
length-scale L, namely, 



3 



D (2) « V'(q) { 1 + 



P 



2 r 



D 



(2) 



V"(qf - 2V'(q)V'"(q) 

VP 



d& « — 

m7 



1 - /3F"(<z)A + /3 2 [lV"{qY + V'(q)V"'(q)] A 2 
1 



1 _ ^-V"(q) 2 \ 2 

4 v ' 



(5) 



The drift and diffusion coefficients specified in (jj) and 
([5]), respectively, provide the systematic extension of the 
leading-order result ([3]) and are applied in the remainder 
to reveal their role in diffusion processes in ratchets. 



III. TRANSPORT IN ADIABATICALLY 
DRIVEN RATCHETS 

We consider a tilted L-periodic potential of the form 

V(q) =V {sin(27rq/i) + a sin [Air{q/L - b)} 
+csin [6ir(q/L — b)]} , 

which lacks an inversion center except for the case a = 
c = 0. An additional random force r\ with zero mean and 
out of equilibrium generates non-vanishing net currents 
and thus directed transport. In calculating these currents 
we pursue the approach followed by Machura et al. in 
[loj and restrict our analysis to dichotomous driving 77 £ 
{—770, +Vo} varying on a time scale much larger than all 
other characteristic time scales of the system (adiabatic 
limit). In this case, the behavior of the system may be 
described by stationary L-periodic solutions -P s t,±r;o (<?) 
of the QSE corresponding to constant currents J± na = 
£Pst,±r) - If ^st.±r;o is normalized to one particle per unit 
cell, this current is for a given bias 77 equal to the mean 
velocity (v) v = LJ n with 



(v)ri 



1 



-0V 



J n dqD 2 (q)-i e-*(«) f Q L+q dye+^y)d(y)^ 



(6) 



Here, the coefficients d and D 2 are given by the respective 
orders in perturbation theory specified in Eqs. ([3]), (U]) 
or ([5]) and the corresponding scaled minimal actions tf> = 
tp/H= Jq dq'Di(q)/D 2 (q) read explicitly 



ft°Hq) = pV{q) -Pm-^r (V'(q) - vY A 



P 2 



(7) 



^\q) = ft°)(q) + ^-V"(q) (V'(q) - V f A 2 , (8) 

^ (2) («)=^ (1) (?). (9) 

Note that the mean velocity © does not depend on D\ 
explicitly. 

The net ratchet-current (v) = [(v) no + (v)_ vo ] /2 is ob- 
tained by expanding (v) n to the respective order in A and 
numerically evaluating the resulting integrals. 
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FIG. 1: (Color online) Net ratchet current (v) vs noise ampli- 
tude r) for the potential |TU} for j3 = 5 and X/L 2 = 0.0025. 
The classical result is represented by a dotted line, the curves 
including leading, locally harmonic and higher-order quantum 
corrections are drawn as short dashed, long dashed and solid 
lines, respectively. 
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FIG. 2: (Color online) Net ratchet current (v) vs noise am- 
plitude 770 for a = 0.4, b = 0.45, c = 0.3, /3 = 2 and 
X/L 2 = 0.0025. Line styles as in fig.Q] 
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FIG. 3: (Color online) Net ratchet current (v) vs inverse tem- 
perature /3o for the same parameters as in fig. [2] and 770 = 5. 
Line styles as in fig. [1] 
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IV. RESULTS AND DISCUSSION 

In order to examine the influence of subleading quan- 
tum fluctuations based on a numerical evaluation of ([5]) , 
it is convenient to introduce dimensionless quantities 
x = q/L, fa = PAV, and W(x) = V(x)/AV where 
the barrier height AV is defined as the difference of the 
maximal and minimal value of the potential for vanishing 
bias. 

We first look at an asymmetric potential with b = c = 
and the following parameter 

W(x) = -0.454 [sin(27ra;) + 0.25 sin(47ra;)] . (f 0) 

In fig.Q]the ratchet current is depicted for various orders 
of perturbation theory together with the classical result. 
Apparently, fluctuations always reduce the current com- 
pared to the classical case. However, the reduction due 
to the leading contributions (of order A) is partially com- 
pensated and the net current pushed back towards its 
classical values when next order terms (local harmonic 
approximation) are considered as well. Only the inclu- 
sion of contributions of order A 2 / L (local anharmonic ap- 
proximation) induces smaller net currents again close to 
the results of the leading fluctuations. The convergence 
of the perturbation series is thus non-monotonous. The 
reason for this is a subtle interplay between higher-order 
contributions in the actions ([7]|9]) and the diffusion coef- 
ficients D2 ■ In local harmonic approximation one always 
has £>2 < and the contribution of fluctuations in 
is reduced by a factor (I — X/3V") compared to that 
in i/j( '. In next order, while ip^ = "ip^, the diffusion 
coefficient is larger than in those regions of 

the ratchet potential, where V" < 0. Accordingly, terms 
produced in order A 2 probe the asymmetry of the ratchet 
weaker than those in order A, and thus soothe quantum 
effects in the net current. In terms of the Euclidian path 
integral for the equilibrium distribution, the local oscilla- 
tor approximation symmetrizes the potential experienced 
by the minimal- action path in a small vicinity around its 



starting and end point q with the tendency to suppress 
quantum phenomena due to asymmetries. 

Let us now analyze this scenario in a more involved 
case a = 0.4, b = 0.45, c = 0.3, V a /AV = 0.372. As 
already reported in [lOj], here leading quantum fluctua- 
tions generate a current reversal compared to the classical 
case, although the absolute value of quantum contribu- 
tions to diffusion and action is small. This reveals the 
sensitivity of the current to even slight variations in the 
ratchet topology. Results for increasing order in pertur- 
bation theory are shown in figs. [2] and [3] Notably, in this 
case data including second order harmonic terms (order 
A 2 ) carry again the typical characteristics of the classi- 
cal net current as a function of the noise amplitude 770. 
Their impact is thus of the same size as that of the leading 
fluctuations. Only when also next order contributions in- 
duced by local asymmetries are taken into account, does 
one obtain results close to leading-order ones. The same 
applies for the situation in fig. [3l where the temperature 
dependence of the current is shown. The conclusion is 
that in the sense of a consistent perturbative treatment, 
it is not justified to neglect formally smaller terms of or- 
der A 2 /L against those of order A 2 . 



V. SUMMARY 

In this paper we have analyzed the impact of sublead- 
ing quantum corrections on net currents in overdamped 
ratchets at very low temperatures. A consistent inclusion 
of these contributions is non-trivial. Beyond the formal 
order of magnitude of the corresponding terms it has to 
take into account also their local symmetry properties in 
position space. Local anharmonicities are essential for 
a convergent perturbative treatment. Along this rule, a 
systematic extension of the QSE provides a powerful tool 
to analyze dynamical properties of strongly condensed 
phase systems at low temperatures. 
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